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Abstract 

In this paper, the characteristic matrix method is employed to theoretically investigate the propagation of 
electromagnetic waves through one-dimensional defective lossy photonic crystals (PCs) composed of negative 
index materials (NIMs) and positive index materials (PIMs). We consider symmetric and asymmetric 
geometric structures with two different types of NIM defect layers at the center of the structure. The effects 
of the polarization and the angle of incidence on the defect modes in the transmission spectra of both 
structures are investigated. The results show that the number of the defect modes within the photonic band 
gap (PBG) depends on the type of the NIM defect layer and is independent of the geometrical structure. 
Moreover, it is shown that the defect mode frequency increases as the angle of incidence increases. This 
property is also independent of the geometry of the structure. The results can lead to designing new types 
of narrowband and multichannel transmission filters. 
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I. INTRODUCTION 



Over the past two and half decades, a new class of material called photonic crystals (PCs) 
has emerged. The PCs originate from theoretical work of Yablonovitch, and John's experimental 
work, which were published almost simultaneously in 1987 .The PCs are artificial dielectric 
or metallic structures in which the refractive index changes periodically in space. This kind of pe- 
riodic structure affects the propagation of electromagnetic waves in the similar way as the periodic 
potential in a semiconductor crystal affects the electron motion by defining allowed and forbidden 
electronic energy bands. Whether or not photons propagate through PC structures depends on 
their frequency. Frequencies that are allowed to travel are known as modes, and groups of allowed 
modes form bands. Disallowed bands of frequencies are called photonic band gaps (PBGs) j], Q]. 
These PBGs are also called the Bragg gaps because they originate from the Bragg scattering in the 
periodic structure. The properties of PBGs in one-dimensional (ID) PCs have been proven to play 
an important role in some potential applications such as photonic devices, optical filters, resonance 
cavities, laser applications, high reflecting omnidirectional mirrors, and the optoelectronic circuits 
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By breaking the periodicity of the conventional PC structure, we will have a defective crystal. 
This can be performed by changing physical parameters, such as changing the thickness of one 
of the layer, adding another medium to the structure, or removing a layer from PCs 9h14||. By 
introducing a layer with different optical properties, localized defect modes, which are also called 
resonant transmission peaks, can be generated within the PBG due to the change of the interference 
behavior of light [lj, [if]] , very similar to the defect states that are generated in the forbidden band 
of doped semiconductors. Upon prediction of the existence of materials with negative refractive 



index (NRI) in 1968 by Veselago [16|], such materials, which have negative permittivity and per- 
meability simultaneously, have received extensive attention for their very unusual electromagnetic 
properties. Negative refractive index materials, or simply negative index materials (NIMs), are 
now also known as metamaterials. Recently, with the possibility of producing metamaterials, PCs 
with metamaterials, called metamaterial photonic crystals (MPCs) have been made. 

In several papers the properties of the defect modes in different ID conventional PC and ID 
MPC structures have been investigated by introdu cing positive or negative indices defects 
351 ] . Following the interesting report by Wu et al. 



231 ] on the properties of the defect modes in 



ID PCs with symmetric and asymmetric geometric structures, here we investigate the properties 
of the defect modes in the transmission spectra of ID defective symmetric and asymmetric lossy 
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MPCs which are composed of negative refractive index material defect layer at the center of the 
crystal. The outline of this paper is as follows. In Section 2, two geometric MPC structures, the 
permittivity and permeability of two types of NIMs, and also the characteristic matrix method and 
its formulation are presented. The numerical results and discussions are given in Section 3, and 
the conclusion is presented in Section 4. 



II. MPC STRUCTURE AND CHARACTERISTIC MATRIX METHOD 



The ID defective MPCs, which are constituted by alternative layers of NIM and positive index 
material (PIM), under study with asymmetric and symmetric structures in air with a defect layer 
at the center of the structures are shown in Figures 1(a) and 1(b), respectively, where NIMs are 
dispersive and dissipative. We assume that layers A and C (defect layer) are NIMs, and layer B 
is a PIM. N is the number of the lattice period, and also di, £j, and m (i = A, B, C) are thickness, 
permittivity and permeability of the layers, respectively. 

The calculations are performed using the characteristic matrix method 3g], which is the most 



effective technique to analyze the transmission properties of PCs. The characteristic matrix for 

(AB) N / 2 C(AB) N / 2 , and symmetric (AB) N / 2 C(BA) N / 2 structures, are given by: 

M[d] = {M a Mb) n/2 M c {M a M b ) n/2 , and M[d] = {M A M B ) N / 2 M C {M B M A ) N / 2 , 

where M A , M B and Mq are the characteristic matrices of layers A, B, and C. The characteristic 



matrix Mj for TE wave is given by 



Mi 



cos7i — sin 7; 



(1) 



-i pi sm7j cos7j 

where, ji = (u>/c) riidiCosOi, c is speed of light in vacuum, 9i is the ray angle inside the layer i 



with refractive index rii, pi = y 'eij ' [ii cos 6i, and cos 9i = J 1 — (?2q sin 2 6q/ii 2 ), in which hq is the 
refractive index of the environment wherein the incidence wave tends to enter the structure. The 
refractive index is given as rii = ±^/Ei[ii [31I, [j?!], where the positive and the negative signs are 
assigned for the PIM and NIM layers, respectively. 

The final characteristic matrix for an N period structure is given by: 

A' 



[M{d)] N = Y[Mi 



(2) 



m n ran 
?7l 2 l m-22 

where m,ij(i,j = 1,2) are the matrix elements of [M(d)] N . The transmission coefficient of the 



i=i 
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multilayer is calculated by: 



* = 7 ^1 V- W 

(mn + mi2 Ps) Po + {m 2 i + m 22 Ps) 



In this equation, pq = uq/ cos 6q and p s = n s j cos8 s , with n s being the refractive index of the 
environment where the wave leaves the crystal with angle 9 S . The transmissivity of the multilayer 
is given by T = (Ps/po)\t\ 2 - The transmissivity of the multilayer for TM waves can be obtained by 



using these previous expressions with p$ = y/ii/ei cos6>j, po = cos Oq/uq, and p s = cos 6 s /n s . 

As mentioned before, our study is based on two different types of NIMs. The permittivity 



and permeability of type-I NIM 
complex, and are defined as [381 . 



layer with negative refracting index in the microwave region are 



£ (f) 1 + Q 9 2 _f2_ i7 f + U 5 2 _ p _ i7 f ( 4 ) 



^ (/) = 1 + 0.9022 _ p _ — f (5) 



where / and 7 are frequency and damping frequency in GHz, respectively. As mentioned in [391 ]. 
for / < 3.13 GHz, the real parts of the permittivity and permeability, e and \i , are simultaneously 
negative (double-negative material). For 3.13 < / < 3.78 GHz, e < 0, but p, > (single-negative 
material extending to the epsilon- negative material), and also for / > 3.78 GHz, both e and /i 
are positive (double-positive material). For type-II NIM layer, we use the Drude model {^(J to 
describe the complex permittivity and permeability, with the results 

M/) = 1.44- 7 ^ y . (7) 

Similarly, / and 7 are respectively frequency and damping frequency, given in GHz. For / < 8.33 
GHz, both e and fi , the real parts of the permittivity and permeability, are negative (double- 
negative material), and for 8.33 < / < 10 GHz, e > 0, but ji < (mu-negative material). 



III. NUMERICAL RESULTS AND DISCUSSION 



Based on the theoretical model described in previous section, the transmission spectrum of the 
lossy defective PC with a NIM defect layer at the center is calculated. The calculations are carried 
out in the region where the real parts of the permittivity and permeability of two types of NIM 

4 
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rs A, and C), e and fx , are simultaneously negative where the zero-n gap will appear 



421 ] . Equations (4) and (5), type-I NIM, are used for the permittivity and permeability 
of layer A. PIM layer (layer B) is assumed to be the vacuum layer with ns = 1. The thickness of 
layers A, B, and C are respectively chosen as (Ia = 6 mm, dg = 12 mm, and dc = 24 mm. Also, 



39] 



the total number of the lattice period is selected to be A = 16 

In the first part, we use equations (4) and (5), type-I NIM, for the permittivity and the per- 
meability of NIM defect layer (layer C) exactly like layer A, and the transmission spectra of two 
different asymmetric and symmetric geometric structures for both polarizations are investigated. 
First of all, the transmission spectra of TE and TM polarized waves for the asymmetric structure 
at various angles of incidence and for 7 = 0.2 x 10 -3 GHz are shown in Figures 2 and 3. As it is 
seen, a single resonant peak appears within the PBG, which corresponds to the type-I NIM defect 
layer. This is in line with the report by Wu et al. 23J] for the defect modes in ID PCs with PIMs 



where only one defect mode for the asymmetric structure has been observed. In addition, the 
figures show that as the angle of incidence increases, the peak height of the defect mode decreases 
for TE waves and increases for TM waves. Moreover, for both polarizations the frequency of the 
defect mode is shifted to the higher frequency as the angle increases. 

In Figures 4 and 5 we have plotted the frequency-dependent transmittance for the symmetric 
structure in TE and TM waves at four different angles of incidence. As it is observed in these figures, 
similar to the asymmetric structure, there is a single defect mode within the band gap again. This 



is in sharp contrast to the work by Wu et al. 



231 ] for the defect modes in ID symmetric PC 



with PIMs where there are two defect modes. It is also worth mentioning that the number of 
defect modes appearing with type-I NIM defect layer does not depend on the geometric structure, 
polarization, and angle of the incidence wave. 

In more detail, as shown in Figure 6, the frequency of defect mode versus the angle of incidence 
for asymmetric (Figure 6(a)) and symmetric (Figure 6(b)) structures have been respectively plot- 
ted. As we mentioned earlier, only one defect mode appears in both asymmetric and symmetric 
structures. Additionally, as seen from the figures, in the asymmetric structure the defect modes 
appear in higher frequencies compared to the symmetric one. Another feature in Figure 6 which 
is worth mentioning is the frequency of the defect mode, which, in TE waves, remains nearly un- 
changed in both geometric structures as the angle of incidence increases. However, in TM waves, 
the frequency increases as the angle of incidence increases. In addition, this increase in symmetric 
structure is higher than that of asymmetric one. 

In the second part, properties of the defect mode of both geometric structures with another type 
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of NIM defect layer at the center (type-II NIM), which was described in Section 2, are investigated. 
The asymmetric and symmetric structures which were used before are modified by replacing type- 
II NIM defect layer (layer C), in which the permittivity and permeability follow the equations (6) 
and (7). The other parameters are kept the same as the first part. 

In Figures 7 and 8, we have respectively plotted the transmission spectra of TE and TM waves 
for the asymmetric structure at various angles of incidence, when the loss factors of two types of 
NIMs, layer A and C, are equal to 7 = 0.2 x 10 -3 GHz. We observe that there are three defect 
modes for both TE and TM waves, as identified by numbers 1, 2, and 3, respectively. Additionally, 
the transmission spectra of TE and TM waves for the symmetric structure have been plotted in 
Figures 9 and 10, respectively. In the symmetric structure, similar to the asymmetric one, there 
are three defect modes again. It can be seen from Figures 7 to 10 that the frequencies of the defect 
modes are shifted to the higher frequency as the incidence angle increases. Besides, these figures 
show that for both asymmetric and symmetric structures, three defect modes appear, which are 
independent of the geometry of the structure. To the contrary, as mentioned by Wu et al. 23], 



in the conventional PCs case for a PIM defect layer, the symmetric geometric structure shows two 
defect modes and the asymmetric one shows one defect mode. 

In more detail, we have plotted the defect mode frequency versus angle of incidence for both 
asymmetric (Figure 11(a)) and symmetric (Figure 11(b)) structures. As seen from these figures, in 
TE waves for both geometric structures, the frequency of defect modes remains nearly unchanged 
with increasing the angle of incidence. On the other hand, frequency of defect modes in TM 
waves increases as the angle of incidence increases. Additionally, similar to the first part, in the 
asymmetric structure the defect modes appear in higher frequencies compared to the symmetric 
one. In contrast to type-I NIM defect layer, frequencies of defect modes, which increase in the 
asymmetric structure, are higher than that of the symmetric structure. Another feature in Figure 
11 is that in the asymmetric structure, defect modes 1 and 2 disappear for angles more than 45° 
and 65°, respectively. Also, defect modes 1 and 2 disappear for angles more than 55° and 70°, 
respectively, in the symmetric structure. To the contrary, Figure 6, by comparing defect mode 3 
in both geometric structures for TM waves, shows that the frequency of the defect mode in the 
asymmetric structure changes more than that of the symmetric one. 

As we observe, in comparison to the first part, which is type-I NIM defect layer, the number of 
defect modes depends on the types of NIMs which are being used for the defect layer. In addition, 
as the figures show, for two types of NIMs, the number of defect modes does not depend on the 



geometric structure. This is in sharp contrast with the work done by Wu et al. 



231 ] where the 
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number of defect modes depends on the asymmetric or symmetric structure. 



IV. CONCLUSION 

This paper has theoretically investigated the properties of defect modes coming from two types 
of NIMs defect layer at the center of ID defective lossy MPCs in symmetric and asymmetric 
geometric structures. Our numerical results show that there exists a single defect mode inside the 
PBG in both asymmetric and symmetric structures for TE and TM polarized waves, when the 
type-I NIM defect layer is considered. On the other hand, the results show that with changing 
the types of NIMs in the defect layer, we observe three defect modes inside the band gap for both 
structures. Positions (frequencies) of the defect modes depend on the polarization of the waves and 
the incidence angles. As the incidence angle increases, the frequency of the defect modes moves 
toward higher frequencies regardless of the geometrical structures and types of NIMs defect layer. 
Additionally, the results of this study bring up the conclusion that the number of defect modes in 
ID asymmetric and symmetric lossy PCs strongly depends on the types of NIM defect. Besides, an 
important result is that by using an NIM in the structure of PCs, the defect modes show different 
behavior from that of the PCs with only PIMs. Finally, it is worthwhile to mention that the 
number of defect modes does not de pen d on the asymmetric or symmetric structures, contrary 
to the work conducted by Wu et al. 



231 ] . Detailed analysis of the defect modes in metamaterial 
photonic crystals with different geometric structures will certainly provide useful information for 
designing new types of narrowband and multichannel transmission filters. 
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FIG. 1: Schematic diagrams of ID defective MPCs embedded in air, with a defect layer, (a) asymmetric 
structure, (b) symmetric structure, where layers A and C are NIMs and layer B is PIM. The number of the 
periods is N and the incidence angle is 0q. 
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FIG. 2: Transmission spectra of TE polarized wave for the asymmetric ID MPC structures with type- 
I NIM defect layer (solid line) and without defect layer (dotted line) for different incidence angles with 
7 = 0.2 x 1(T 3 GHz. 
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FIG. 3: Transmission spectra of TM polarized wave for the asymmetric ID MPC structures with typc- 
I NIM defect layer (solid line) and without defect layer (dotted line) for different incidence angles with 
7 = 0.2 x 1(T 3 GHz. 
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FIG. 4: Transmission spectra of TE polarized wave for the symmetric ID MPC structures with type-I NIM 
defect layer (solid line) and without defect layer (dotted line) for different incidence angles with 7 = 0.2 x 1CP 3 
GHz. 
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FIG. 5: Transmission spectra of TM polarized wave for the symmetric ID MPC structures with type-I NIM 
defect layer (solid line) and without defect layer (dotted line) for different incidence angles with 7 = 0.2 x 10~ 3 
GHz. 
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FIG. 6: Frequency of the defect modes in (a) the asymmetric and (b) the symmetric ID PMCs as a function 
of angle of incidence for both polarizations, with 7 = 0.2 x 10~ 3 GHz. 
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FIG. 7: Transmission spectra of TE polarized wave for the asymmetric ID MPC structures with typc- 
II NIM defect layer (solid line) and without defect layer (dotted line) for different incidence angles with 
7 = 0.2 x 1(T 3 GHz. 



o 

I 
a 



I 

B 

i 




1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



1.5 2.0 2.5 

Frequency-[GHz] 



, = 25° 

















2 


3 




9 a = 50° 




u 


y 





2.0 

Frequency- 




FIG. 8: Transmission spectra of TM polarized wave for the asymmetric ID MPC structures with type- 
II NIM defect layer (solid line) and without defect layer (dotted line) for different incidence angles with 
7 = 0.2 x 10~ 3 GHz. 
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FIG. 9: Transmission spectra of TE polarized wave for the symmetric ID MPC structures with typc-II NIM 
defect layer (solid line) and without defect layer (dotted line) for different incidence angles with 7 = 0.2 x 10~ 3 
GHz. 
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FIG. 10: Transmission spectra of TM polarized wave for the symmetric ID MPC structures with typc- 
II NIM defect layer (solid line) and without defect layer (dotted line) for different incidence angles with 
7 = 0.2 x 10~ 3 GHz. 
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FIG. 11: Frequency of the defect modes in (a) the asymmetric and (b) the symmetric ID PMCs as a function 
of angle of incidence for both polarizations, with 7 = 0.2 x 10~ 3 GHz. 
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